The effect of Landau-Zener dynamics on phonon lasing 
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Optomechanical systems couple light to the motion of nanomechanical objects. Intriguing new 
effects are observed in recent experiments that involve the dynamics of more than one optical mode. 
There, mechanical motion can stimulate strongly driven multi-mode photon dynamics that acts 
back on the mechanics via radiation forces. We show that even for two optical modes Landau- 
Zener- St ueckelberg oscillations of the light field drastically change the nonlinear attractor diagram 
of the resulting phonon lasing oscillations. Our findings illustrate the generic effects of Landau-Zener 
physics on back-action induced self-oscillations. 
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The exploration of nanomechanical objects and their 
interaction with light constitutes the rapidly evolving 
field of optomechanics (see [1 for a recent review). The 
key element of any optomechanical system is a laser- 
driven optical mode whose resonance frequency shifts in 
response to the displacement of a mechanical object. The 
photon dynamics conversely acts back on the mechanics 
in terms of a radiation pressure force. These dynamical 
back-action effects, mediated by the light field, can cool 
or amplify mechanical motion, and even drive the system 
into a regime of self-induced mechanical oscillations [2]-[8] 
akin to lasing. 

An exciting new development has introduced optome- 
chanical setups with multiple coupled optical and vibra- 
tional modes [9HTT]. These systems allow one to real- 
ize sophisticated measurement schemes p!2HT4j , to study 
collective phenomena p^, or to mechanically drive co- 
herent photon dynamics [16]. For applications, they fur- 
thermore stimulate prospects of integrated optomechan- 
ical circuits [Fr, TS] . Recently, phonon lasing for an op- 
tomechanical setup involving a tunable optical two-level 
system has been demonstrated [19]. Hence, implement- 
ing a nanomechanical analog of a laser has finally been 
achieved [20 . Here, we develop the fully nonlinear theory 
of phonon lasing (self-induced mechanical oscillations) in 
such multimode optomechanical setups. In particular, 
we will point out that the mechanical oscillations may 
induce Landau-Zener physics with respect to the optical 
two-level system, and that this has a strong effect on the 
dynamical back-action. The resulting phenomena drasti- 
cally change the nonlinear attractor diagram, i.e. the re- 
lation between the mechanical lasing amplitudes and the 
experimentally tunable parameters. We will refer to an 
existing optomechanical setup [9] where our predictions 
can be verified experimentally. However, most of our 
analysis and discussion are applicable to the quite generic 
situation where self-induced oscillations are pumped by 
a parametrically coupled, driven two-level system. Our 
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FIG. 1. (a) Setup. A moveable membrane, placed inside a 
cavity, couples two optical modes a^, an via transmission, 
(b) Optical resonance frequencies vs. displacement. The 
membrane's displacement linearly changes the bare modes' 
frequencies (dashed). Due to the photon coupling g there is 
an avoided crossing for the resonance frequencies uj± (black). 
Mechanical oscillations x{t) = Acos(Qt) + Xa periodically 
sweep the system along the photon branches (red), (c) Cav- 
ity resonance frequency (jJ±{x{t)) depending on time. For 
non-adiabatic sweeps through the anti-crossing, repeated LZ 
transitions (highlighted regions) split the photon state. Af- 
ter each passage, the two contributions gather a phase differ- 
ence that leads to subsequent interference. The resulting LZS 
oscillations in the light field act back on the mechanics via 
the radiation pressure force, (d) For sufficiently large back- 
action-induced anti-damping, the system enters a regime of 
mechanical self-oscillations (phonon lasing). 



findings thus are also relevant for nanomechanical struc- 
tures or microwave modes whose oscillations are ampli- 
fied by coupling to, e.g., current-driven double quantum 
dot setups, superconducting single-electron transistors, 
or Cooper-pair boxes [21]. 

We consider the system depicted in Fig. [l^. A dielec- 
tric membrane is placed in the middle between two high- 
finesse mirrors [9]. Transmission through the membrane 
couples the optical modes of the left and right half of the 



cavity, respectively. Focussing on two nearly degenerate 
modes, the Hamiltonian of the cavity reads 



(1) 



Here, Oj^cll {clj^clr) is the photon number operator of the 
left (right) cavity mode, ujq is the modes' frequency for 
X = (where the two modes are degenerate), and 21 is the 
length of the full cavity. The membrane's displacement 
X linearly changes the modes' bare frequencies, while the 
optical coupling g leads to an avoided crossing for the 
system's two optical resonances, uj± = \/ ^ {ujqx/IY 
[Fig. Thus, mechanical oscillations x{t) periodically 
sweep the system along the hyperbola branches uj±. 

We focus on the experimentally accessible, non- 
adiabatic regime [12] [16] where fast periodic sweep- 
ing through the avoided crossing results in consecutive 
Landau-Zener (LZ) transitions [221123]. For a photon in- 
serted into the left mode, the first transition splits the 
photon state into a coherent superposition, the two con- 
tributions gather different phases and interfere the next 
time the system traverses the avoided crossing [Fig. [l]^]. 
For a two-state system, the resulting interference pat- 
terns are known as Landau- Zener-Stueckelberg (LZS) os- 
cillations [2l]. These have been demonstrated in many 
setups, ranging from atomic systems |25h22| to quan- 
tum dots and superconducting qubits [28H3T] . In all of 
these situations, LZS effects are produced by a fixed ex- 
ternal periodic driving. In contrast, here we address 
the case where LZS oscillations act back on the mech- 
anism that drives them (i.e. the mechanical motion), via 
the radiation pressure force. We will see that LZS in- 
terference strongly influences this back-action force and 
thereby drastically affects the mechanical self-oscillations 
that occur when this force overcomes the internal friction 
[Fig. [l]i] . More generally, the following discussion thus 
illustrates the effect of LZS dynamics on back-action in- 
duced instabilities. 

Given the radiation pressure force Frad — —dHcav/dx^ 
the coupled equations of motion for the displacement x{t) 
and ai{t) {i = L^R), read 

X = Ao{d'^^aL - o.'^rClr) - n'^{x - xo) - Fx + i{t), (2) 
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where we used input-output theory for the light fields and 
set ^0 = J^o/lm. The membrane has a mechanical fre- 
quency 1], an intrinsic damping rate F and a rest position 
xq. Photons decay at a rate out of the cavity. We as- 
sume the left mode to be driven by a laser at frequency 
ujl] the input fields bl^{t) contain this contribution. In 
the following, we will consider purely classical (large- 
amplitude) nonlinear dynamics and replace the operators 



(a) 













it 


< 

X : ^/'^ 






(b) . 


^ _L ^ 


(c) 




▼ 






-1.0 -0.5 0.0 0.5 1.0 

Xal^ 



OlR . 



FIG. 2. (a) Effective optomechanical damping Topt for given 
mechanical oscillations x{t) = Acos(Qt) -\- Xa as a function 
of mean position Xa and laser detuning A^- Parameters: 
A/Q = 0.5, g/Q = 0.2, k/Q = 0.1. Mechanical sidebands 
(dashed), displaced by multiples of Q, show cooling (blue; 
Fopt > 0) and heating (red; Fopt < 0). iFoptl is largest if 
the optical modes' frequency difference is in resonance with 
the mechanical frequency Q; position (b) and (c). For finite 
amplitude, this yields an Autler-Townes splitting (see circled 
regions). Topt in units of 2ujoPin/mQ^l^ . (b) Creation (heat- 
ing) or (c) destruction (cooling) of a phonon upon transferring 
a photon from left to right, (d) At the degeneracy point, the 
bare optical frequencies are swept past each other in an oscil- 
latory fashion (cf. Fig. [ij)). 



ai{t) by the coherent light amplitudes ai{t). The classical 
input fields then read = 0, = e-'"^^^ ^fP~JnLjl 
where Pin is the laser input power, and the mechanical 
Langevin force will be neglected (^ ~ 0). For conve- 
nience, we define the laser detuning /S.l = uol — uoq. 

The radiation pressure force gives rise to a time- 
averaged net mechanical power input (Fradx). In anal- 
ogv to the intrinsic friction F, see Eq. Il), we can define 
(Fradx) = —mV opt{x^) such that we obtain an effective 
optomechanical damping rate 

^0 



opt 
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For Topt > {Topt < 0) the light-field interaction 
damps (anti-damps) the mechanics. For given oscilla- 
tions x{t) = Acos{Qt)-\-Xa, Topt can be calculated via the 
periodic light field dynamics aL{t), OLR{t) that is found 
by solving Eq. (|3|; see also Eq. (|6| further below. Note 
that our Fopt is amplitude-dependent, and the usual lin- 
earized case [l3l[32] is recovered for A ^ 0. In the follow- 
ing we will express displacement in terms of frequency, 
x{t) = {uJo/l)x{t) (see Eq. ([l])); likewise for A, Xa- 

Fig. [2^ shows results for Topt in this setup, at moder- 
ate amplitudes A. Optomechanical damping and heating 
is largest if the optical modes' frequency difference is in 
resonance with the mechanical frequency [8] [19]. In 
this case, photon transfer from the laser-driven left mode 
into the right one involves absorption (or emission) of a 
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phonon, that yields strong mechanical heating (or cool- 
ing), see Fig. [2]3-c. For finite amplitudes, we observe an 
Autler-Townes (AT) splitting [33 that scales as 2gA/rt 
[El. Given 

^opt^ we now turn to discuss back-action 
driven mechanical self-oscillations (phonon lasing) of the 
membrane. 

For suitable laser input powers, the radiation pressure 
force only weakly affects the mechanics over one oscilla- 
tion period and the mechanics approximately performs 
sinusoidal oscillations at its unperturbed eigenfrequency 
x{t) = Aco^{Vtt) + Xa- The possible attractors of 
the dynamics (A, x^) have to meet two conditions [5[|6]. 
First, the time- averaged total force must vanish: {x) =0. 
Second, the overall mechanical power input due to radi- 
ation pressure must equal the power loss due to friction, 
{xx) = 0. From Eq. [2j the power balance {xx) = is 
equivalent to 



^opt{^-) ^a) — r. 



(5) 



The force balance {x) = yields (Frod(^)) = 
mft'^ {xa — xo), i.e. the radiation pressure force displaces 
the membrane's average position Xa from its rest posi- 
tion xq. In general, one solves the force balance to find 
= Xa{A,xo) and uses this to calculate ToptiA^Xa) 
[5l|6]. For high quality mechanics {^/T 1), the power 
balance (Eq. ([5|) is met for weak radiation pressure forces 
where Xa — xq. For clarity, we will focus on this case. 
Otherwise, attract or diagrams get deformed slightly [5]. 

Fig. [3^ displays the effective optomechanical damping 
Topt depending on laser-detuning A^, and amplitude A. 
The structure of this diagram is drastically different from 
the standard case with one optical mode [5j [6] . There are 
"ridges" of high Fopt which display an oscillatory shape 
(clarified in the inset). 

A physical understanding of Fig. |3] can be found from 
the general structure of the light field dynamics that en- 
ters the optomechanical damping, Eq. Q. For given me- 
chanical oscillations x{t) = Acos{Qt) + x^, the formal 
solution to Eq. Q can be expressed as 



f 



Gi{t,t')e-^'-*-'">/^e-'^^*'dt' 



(6) 

where the Green's function Gi(t^t') describes the ampli- 
tude for a photon entering the left mode at time t' and to 
be found in the left or right one {i = L^R) time t. From 
Eq. ^ Gi(t,t')\s found to be Gi{t,t') = ai{t,t')e-''t'^^'^ 
where (j){t') = (A/ft) sm{rtt^) and ai{t^t') is a solution to 
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with t>t' and initial condition aR{t' ^ t') = 0, aL{t\ t') = 
1. Thus, the internal photon dynamics between the two 
modes ai(t^t') is expressed in terms of a two- level sys- 
tem with a time-dependent coupling ge^'^^^'^\ With = 




FIG. 3. Attractor diagram for phonon lasing oscillations 
(regime Q > 2g). (a) Effective optomechanical damping Vopt 
as a function of laser detuning A l and oscillation amplitude 
A, for a membrane positioned at the degeneracy point Xa = 0; 
other parameters as in Fig. [2] For large amplitudes, the in- 
terference between consecutive LZ transitions (Fig. [l]:) leads 
to LZS oscillations. They result in ridges of high Fopt, whose 
oscillatory shape can be understood via the Floquet eigen- 
values e± (middle panel) for the periodic light field dynam- 
ics. The ridges are located at Ai, = mQ + tj{A)^ where 
e±{A) ^ zb^Jo(2^/Q) involves a Bessel function. Dashed 
lines indicate |Al| = A. (b) Blow-up of framed region in 
(a). The contour lines at Vopt{A^Xa) — — F (Eq. [5]) denote 
possible attractors (allowed amplitude values: solid - stable 
/ dashed - unstable) for the mechanical self-oscillations gen- 
erated by back-action; plotted for two different values of — F, 
as indicated. (F, Topt in units of 2uQPin / mOi^l^) 



{clr^cll)^^ Eq. ^ is the Schrodinger equation including 
a time-periodic Hamiltonian, H{t + T) = H(T). In this 
case it is appropriate to consider the time-evolution op- 
erator for one period, ilj{t' + T) = U{T)tlj{t')^ and its 
two eigenvalues, the so-called Floquet eigenvalues e±: 
U{T)x± = exp(-ze±T)x±. U{T) is obtained by inte- 
grating Eq. 

Using Floquet theory [3T, we find the gen- 
eral structure of the Green's function Gi{t^t') = 
E,,n,n'^r'"''e"'^^"'""''^e-^^^(^-^'), where Cf'"''^' are 
time-independent coefficients. Then, via Eq. (|6| we ob- 
tain pronounced resonances in Fopt located at Ai, = 
mVt + e±{A)^ corresponding to the ridges in Fig. (3] The 
interference between consecutive LZ transitions renor- 
malizes the coupling between modes in terms of Bessel 
functions J^: ge'^''^^^^ = g^^ Jn{2A/n)e'''^^ (Eq. [t]). 
This results in an oscillatory modulation of the Floquet 
eigenvalues e±{A). At certain amplitudes, these vanish 
due to total destructive interference, see Fig. |3^. The 
oscillatory shape of the ridges in Fopt then directly de- 
termines the attractor diagram for the self-induced oscil- 
lations, via the power balance equation ([5|, see Fig. [sJd. 

Regarding the global structure of Fig.|3^, Fopt tends to 
be large near A^ = ±A (dashed lines). This is because 
then the left mode gets into resonance with the laser at 
the motion's turning point. For larger amplitudes, we 
recover the predictions for the standard optomechanical 
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FIG. 4. (a) Overview of the parameter space. Multimode 
dynamics leads to effects beyond the standard scenario when 
the optical splitting can be resolved: 2g > k, (colored region). 
Parameter set (1) corresponds to the one in Figs. [2] [s] set (2) 
is considered in (b-c). (b) Effective optomechanical damping 
Topt for given mechanical oscillations x(t) = Acos{Qt)-\-Xa, as 
a function of mean position Xa and laser detuning Al (com- 
pare Fig. [2]). Parameters: A/Q = 1.5, g/Q = 2.3, k/Q = 0.2. 
(c) Attractor diagram. Effective optomechanical damping 
-L opt as a function of laser detuning A l and oscillation ampli- 
tude A for a membrane positioned at the degeneracy point, 
Xa = 0. Further parameters as in (b). The solid contour 
line FoptiA^Xa) = — F indicates the stable attractors for self- 
induced oscillations. Green (thick) lines show the asymptotic 
behavior. (Topt in units of 2iJoPin/mQ^l^) 



setup [5 (checkerboard in Fig. |3^). 

So far, we discussed dynamical back-action effects for 
parameters where the mechanical frequency is larger than 
the optical splitting, Q > 2g [Fig. 2|3 



In general, the 

parameter space can be subdivided as shown in Fig. [4^. 
Multimode dynamics that goes beyond the standard sce- 
nario ^ i6j can only be observed if the photon lifetime 
inside the cavity is larger than the timescale for pho- 
tons to tunnel between modes, 2g > n (colored region. 
Fig. [4^). Otherwise, photons inserted into the left mode 
decay before the second mode affects the dynamics and 
we recover the standard results [S] |6]. Within the new 
region (colored in Fig. [4^), the most interesting regime is 
where mechanical sidebands can in fact be resolved, i.e. 
K < Vt. Above, we had focussed on the sector 2g < Vt 
within this regime. Now Fig. [4]3 displays Topt in the op- 
posite sector where 2g > Vt. Here, several mechanical 
sidebands lie within the avoided crossing. With respect 
to self-induced mechanical oscillations, these sidebands 
and their interaction yield an intricate web of multistable 
attractors, see Fig. The global asymptotics of these 
structures (green lines) can be found from the quasistatic 
approximation, i.e. from the tim e- average d transition 
frequency: Al = 2(cj+(t)) = 4.^ g'^ ^ A^E{tt /2,k) /it , 



To conclude, we have investigated self-induced me- 
chanical oscillations (phonon lasing) in a multimode op- 
tomechanical system. The mechanical motion drives 
Stueckelberg oscillations in the light field of two coupled 
optical modes, and this drastically modifies the attrac- 
tor diagram. The additional influence of quantum (and 
thermal) noise could be analyzed along the lines of [5] [6] . 
Our example, which can be realized in present optome- 
chanical setups, illustrates the potential of Landau-Zener 
physics to appreciably alter lasing behavior. 
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where k = \j / (^^ + £^(|, /c) is the complete 

elliptic integral of the second kind. 
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